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Abstract 



A search for large extra spatial dimensions via virtual-graviton exchange in the dipho- 
ton channel has been carried out with the CMS detector at the LHC. No excess of 
events above the standard model expectations is found using a data sample collected 
in proton-proton collisions at = 7 TeV and corresponding to an integrated lu- 
minosity of 36 pb^^. New lower limits on the effective Planck scale in the range of 
1.6-2.3 TeV at the 95% confidence level are set, providing the most restrictive bounds 
to date on models with more than two large extra dimensions. 
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1 Introduction 

Compact large extra dimensior\s (ED) are an intriguing proposed solution to the hierarchy 
problem of the standard model (SM), which refers to the puzzling fact that the fundamental 
scale of gravity Mpi ~ 10^^ GeV is so much higher than the electroweak symmetry breaking 
scale ~ 10'^ GeV. With such a difference in scales, it is difficult to protect the Higgs boson mass 
from radiative corrections without a very high degree of fine-tuning. The original proposal 
to use ED to solve the hierarchy problem was presented by Arkani-Hamed, Dimopoulos, and 
Dvali (ADD) flj |2|. They posited a scenario wherein the SM is constrained to the common 
3+1 space-time dimensions, while gravity is free to propagate through the entire multidimen- 
sional space. Thus, the gravitational flux in 3+1 dimensions is effectively diluted by virtue of 
the multidimensional Gauss's Law. The fundamental Planck scale Md is therefore related to 
the apparent scale Mpi according to the formula M'^°^^ = Mp^/ r"^o^ where r and ngo are the 
size and number of the extra dimensions, respectively. A review of current limits on the ADD 
model can be found elsewhere 

Phenomenologically, this large ED scenario results in s-channel production of massive, virtual 
Kaluza-Klein (KK) graviton states, which can decay into two photons ||4|. Decays to fermions 
are suppressed relative to photons because the graviton is spin-2 and the fermions cannot be 
produced in the s wave. Because the KK gravitons propagate through the compact extra dimen- 
sions, their wavefunction must satisfy periodic boundary conditions, which result in discrete 
energy levels with modal spacing on the order of the inverse ED size, from 1 meV to 100 MeV. 
This results in an apparent continuum spectrum of diphotons, rather than distinct resonances. 

Summing over KK modes results in a divergence in the cross section, so an ultraviolet (UV) 
cutoff scale Ms is needed. This scale is related to — ^but potentially different from — the fun- 
damental Planck scale Md, as the precise relationship depends on the UV completion of the 
effective theory. The effects of virtual-graviton production on the cross section are parameter- 
ized by the single variable rjc = T I Mg, where T is an order-unity dimensionless parameter, 
for which several conventional assumptions exist: 



T = 1 (Guidice, Rattazzi, and Wells, GRW [5]), (1) 
log(^) ifnED = 2 




if Med > 2 



(Han, Lykken, and Zhang, HLZ lH), (2) 

(3) 

where is the center-of-mass energy of the hard parton-parton collision. We note that the 
HLZ convention contains an explicit dependence on the number of extra dimensions. 

Searches for extra dimensions via virtual-graviton effects have been conducted at HERA, LEP, 
and the Tevatron (Refs. O IH contain recent reviews of these searches). The most stringent 
previously published limits on Ms come from the DO measurements in the dijet |9| and dipho- 
ton plus dielectron [lOJ channels, which exclude values of Ms lower than 1.3—2.1 TeV at 95% 
confidence level (CL), depending on nEo- 

In this paper, we present a search for virtual-graviton contributions in the diphoton final state, 
using a data sample corresponding to an integrated luminosity of 36 pb^^, collected in pp 
collisions at a/s = 7 TeV at the CERN Large Hadron Collider (LHC) with the Compact Muon 
Solenoid (CMS) detector. 
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3 Event Reconstruction and Selection 



2 The CMS Detector 

CMS is a general-purpose detector designed to study proton collisions at the LHC and is de- 
scribed in detail in Ref. IfTT) . The detector consists of an all-silicon tracker, an electromagnetic 
calorimeter (ECAL), and a hadronic sampling calorimeter (HCAL), all contained inside a large- 
bore 3.8 T superconducting solenoid. In the central region, the tracker consists of three layers 
of silicon pixel detectors, followed by 10 layers of single- and double-sided silicon-strip detec- 
tors. The calorimeter towers are projective and finely segmented, with A/y ^ A<p ^ 0.087 in the 
central region. Moreover, each tower consists of a five-by-five transverse grid of ECAL crystals 
(At] 0.0174), allowing precise reconstruction of the e/ 7 position and energy. Here, the 

pseudorapidity f/ is defined as — In (tan j), where 6 is the polar angle with respect to the di- 
rection of the counterclockwise beam, and (p is the azimuthal angle. Beyond the solenoid there 
are four layers of muon detectors, which are interspersed throughout the steel return yoke of 
the magnet. The instantaneous luminosity is measured with a relative uncertainty of 4% using 
information from forward hadronic calorimeters [12J. 

The CMS trigger system consists of two levels. The first level (LI), composed of custom hard- 
ware, uses irvformation from the calorimeters and muon detectors to select the most interesting 
events for more refined selection and analysis at a rate of up to 80 kHz. The software-based 
High Level Trigger further decreases the rate to a maximum of ~ 300 Hz for data storage. 
Events for the control samples used in this analysis were collected through a single-photon 
trigger, where the photon was required to have a transverse energy Ej = E sin 9 of at least 
20 or 30 GeV, depending on the data collection period. Events in the signal sample were also 
initially collected through a single Ej > 20 GeV photon trigger. However, as the instantaneous 
luminosity increased, we switched to a double-photon trigger requiring both photons to have 
Et > 22 GeV. 

3 Event Reconstruction and Selection 

We first require that an event be consistent with a pp collision and have at least one well- 
reconstructed primary vertex ||T3| . We then reconstruct photons with Ej > 30 GeV in the 
ECAL barrel fiducial region {\rj\ < 1.44) by clustering electromagnetic energy depositions in 
the ECAL. The ECAL clusters are five crystals wide in t] and a variable length in (p to capture 
associated electromagnetic energy from possible photon conversions in the tracker. If hits are 
present in the pixel detector consistent with an electron track whose momentum and location 
is similar to the energy and location of the ECAL cluster, then the cluster is rejected as a photon 
candidate. In 2010 collision data, the ECAL has an energy resolution better than 1% in the 
barrel for unconverted photons with Ej > 20 GeV IImH . 

Hadronic jets can be misidentified as photons when their leading hadron is a hard or rj. 
We reduce the misidentification rate by placing the following restrictions on the isolation of 
the cluster: (i) the hadronic energy within AR < 0.15 of the cluster must be less than 5% of 
its electromagnetic energy; (ii) the scalar sum of the transverse momentum of tracks, Epx/ 
associated with the primary event vertex surrounding the cluster within a hollow cone of 
0.04 < AR < 0.40 must be less than 2.0 GeV + O.OOlEj/ where Ej is the photon transverse 
energy (a rectangular strip of At] x A^ = 0.015 x 0.400 is excluded from the track pj summa- 
tion to allow for photons that convert into e^e^ pairs); (iii) the EEj of ECAL energy surround- 
ing the cluster within 0.06 < AR < 0.40 (and excluding a strip of Arj x Acp = 0.04 x 0.400) 
must be less than 4.2 GeV + 0.006Et; and (iv) the EEj of HCAL energy surrounding the cluster 
within 0.15 < AR < 0.40 must be less than 2.2 GeV + 0.0025Et. Here AR = y^A(p^ + A//2 = 
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\/{<p — <P-yY + ('/ ~ ViY- ^^^o require that the shower shape in rj, £7^,^, be consistent with a 
photon. The c^y,; variable is a modified second moment of the electromagnetic energy cluster 
about its mean t] position, defined in Ref. [15J. Topological and timing criteria suppress noise 
present in the ECAL [14, 16 1. 

We reconstruct two photons using the selection described above and require that the invariant 
mass of the two photons satisfies M-y-y > 60 GeV. The invariant mass and photon pseudora- 
pidity selection criteria are optimized to produce the highest sensitivity for values of Ms and 
Med to which the present data is sensitive. Alternating between fixing the M^y-y requirement 
and floating the \rj\ requirement — and vice versa — results in a final choice of |//| < 1.44 and 
M-y-y > 500 GeV. The optimal choice for | rj \ was very close to the ECAL barrel-endcap bound- 
ary, so we chose the edge of the transition region into the boundary for simplicity. This defines 
the signal region. The intervals 60 < M^-y < 200 GeV and 200 < M^^y < 500 GeV define the 
control and intermediate regions of the data, respectively. 

The photon reconstruction and identification efficiency is measured in Monte Carlo (MC) sim- 
ulation and corrected using a data/MC scale factor of 1.010 ± 0.012 derived from studying 
Z — > e+e^ events. The final efficiency is roughly constant as a function of the photon Ej arid rj. 
The efficiency for an £x > 30 GeV photon with |?7 1 < 1.44 is (87.8 ± 2.3)%, where the dominant 
systematic uncertainty is chosen to cover the variation as a function of and //. Therefore, the 
corresponding diphoton reconstruction and identification efficiency is (77.1 it 4.5)%. 

4 Signal and Background Estimation 

We simulate ED in the ADD model using the SHERPA (vl.1.2) f\7\ MC generator, which sam- 
ples different operating points in Ms and n^D, followed by a fast parametric simulation of the 
CMS detector IflS) . A fast simulation is adequate for describing multiphoton final states and 
has been extensively validated using full simulation of the detector via GEANT4 [19J. The simu- 
lation includes both SM diphoton production and signal diphoton production via KK-graviton 
exchange, in order to account for the interference effects. We use CETQ6l1 [20J parton distri- 
bution functions (PDF) in the simulation. The leading order (LO) SHERPA cross sections are 
multiplied by a next-to-leading order (NLO) K factor of 1.3 ± 0.1 [21J, where the systematic 
uncertainty covers the variation of the K factor with the diphoton mass. Additionally, a 1.5% 
relative uncertainty on the signal acceptance is included to account for uncertainty due to the 
PDF. 

Backgrounds due to the mimicking of a photon signal by a jet are small in the signal region. 
There are two sources of these backgrounds from isolated photon misidentification which we 
consider: multijet production and prompt photon production (i.e., photons from 7 + jets). In 
particular, we measure a misidentification rate, defined as the ratio of the number of isolated 
photons to non-isolated photons in a sample, where the non-isolated photons are selected sim- 
ilarly to the isolated photons except that they fail one of the isolation or shower-shape criteria. 
The definition of the numerator and denominator are given so that the two sets are exclusive 
to each other. The misidentification rate is measured in an EM enriched sample and is then ap- 
plied to the observed events with one or more non-isolated photons, resulting in a prediction 
of the dijet and 7 + jet backgrounds. 

Because a given control sample in which we measure the misidentification rate may contain 
some number of real, isolated photons that "contaminate" the misidentification-rate numera- 
tor, we correct for the numerator purity on a sample-by-sample basis. This is done by releasing 
the arm requirement and fitting the numerator sample for the fraction of prompt photons us- 
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4 Signal and Background Estimation 



ing one-dimensional probability density histograms ("templates") in o^^fj. The signal template 
is constructed from MC simulation, and the background template is constructed by using in- 
verted isolation criteria. The measured misidentification rate falls from 28% at Ej = 30 GeV to 
2% at Ej = 120 GeV. We use two other complementary techniques (using converted photons 
and an isolation template to estimate prompt-photon contamination) to bound the misiden- 
tification rate and apply a conservative 20% systematic uncertainty which is the dominant 
uncertainty on the background estimation. 

The diphoton background is computed with the SHERPA MC program and then rescaled by an 
NLO K factor of 1.3 ETHm. This K factor is alternatively derived with DIPHOX |23|, wherein 
we observe that the K factor decreases slowly as a function of diphoton invariant mass and 
stabilizes in the range of interest. We therefore use a K factor of 1.3 ± 0.3 for the diphoton 
background to cover its variation as a function of My^ throughout the control, intermediate, 
and signal regions (M-y-y > 60 GeV). We observe (expect) 440 (374 ± 51) events with Mjj > 
60 GeV and zero (0.30 ± 0.07) with Myj > 500 GeV. 
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Figure 1: Distributions in t] and £7- for the leading and sub-leading photons. Points with error 
bars represent observed data; the solid histogram corresponds to the expected background 
from control regions. Shading corresponds to the systematic uncertainty on the background 
expectation. 
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Figure [T] shows the pseudorapidity and transverse energy distributions for the photons with 
the largest and second-largest transverse energy in an event, along with the background pre- 
dictions. Figure |2] displays the invariant mass distribution for each of the backgrounds as well 
as the observed data, with the optimized tj requirement. Beyond M-y-y > 200 GeV, we see 12 
events, consistent with 9.0 ± 1.5 expected background events. 




100 200 300 400 500 600 700 800 900 1000 

M,, (GeV) 

Figure 2: Observed data (points with error bars) and background expectations (filled solid 
histograms) as a function of the diphoton invariant mass. Photons are required to be isolated, 
with £x > 30 GeV and \y]\ < 1.44. Also shown with dashed lines are the signal distributions for 
two sets of model parameters. Shaded bands around the background estimation correspond to 
systematic uncertainties. The last bin is an overflow, including the sum of all contributions for 
M^^ > 1.0 TeV. 

Table [l] presents the number of observed data events and expected number of background 
events in the control, intermediate, and signal regions, respectively. This table corresponds di- 
rectly to the plot in Fig. |2] In the control region, we find that the data are consistent with the 
background expectation within the systematic uncertainty. The systematic uncertainty on the 
the total background takes into account the correlations between dijet and 7 + jet backgrounds 
arising because both depend on the same misidentification rate. The relative combined uncer- 
tainty on the backgrounds in the signal region is 22%, due nearly entirely to the diphoton NLO 
X-factor uncertainty. 
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5 Results 



Table 1: Data yields and background expectations for reconstructed diphoton invariant mass 
ranges. Full systematic uncertainties have been included. 



Process 


60 < Myy < 200 GeV 


200 < Mjy < 500 GeV 


500 GeV< Myy 


Dijets 


70 ±28 


0.5 ±0.2 


0.0009 ± 0.0004 


7 + Jets 


145 ±7 


2.3 ±0.3 


0.016 ± 0.003 


Diphotons 


150 ± 35 


6.2 ±1.4 


0.29 ± 0.07 


Total Backgrounds 


365 ± 49 


9.0 ±1.5 


0.30 ± 0.07 


Observed 


428 


12 






Table 2: Summary of systematic uncertainties. 





Central Value 


Relative Uncertainty 


Luminosity 


36 pb-^ 


4.0% 


Background diphoton K factor 


0.30 Events 


23% 


Signal Efficiency 


77.1% 


6.0% 


Signal diphoton K factor 


1.3 


7.7% 



5 Results 

We perform a counting experiment in the signal region {Mjy > 500 GeV) and set 95% CL upper 
limits on the quantity 

S = (t^total - CT-Sm) X^xA, (4) 

where c^total represents the total diphoton production cross section (including both signal, SM, 
and interference effects), and usm represents the SM diphoton production cross section, only. 
We indicate the signal branching fraction to diphotons by /3 and the signal acceptance by A. We 
utilize a standard Bayesian approach ||3H24| with a flat prior chosen for the signal cross section 
and log-normal priors for the nuisance parameters (integrated luminosity, signal efficiency, 
and background). The likelihood is constructed from the Poisson probability to observe N 
events, given S, the signal efficiency of (77.1 ± 4.5)%, the expected number of background 
events (0.30 ± 0.07), and the integrated luminosity C = (36 ± 1) pb^^ [12J. Table [2] summarizes 
the systematic uncertainties used as inputs to the limit calculation. 

The observed (expected) 95% CL upper limit on S is 0.11 (0.13) pb. We translate the limit on S 
into a limit on the parameters of the ADD model using the following technique. Because the 
effects of virtual-graviton exchange interfere with the SM diphoton production, generally, we 
expect the total diphoton cross section to have the following form: 

C^total = C^SM ± f]G C^int ± Vc '^ED/ (5) 

where c%t arid (7ed refer to the contributions to the total cross section from interference and di- 
rect ED effects, respectively. Here, the dimensionful }]g parameter specifies the strength of the 
ED effects and is related to Ms through Eqs. ([l]-Q. Consequently, we parameterize S as a func- 
tion in the parameter j/g- For the HLZ = 2 case, rjc depends on the signal spectrum due to 
an explicit s dependence in Eq. (|2]). Therefore, in this case we parameterize S as a function of 
1/Mg and translate the limit on this parameter into a limit on Mg. 

The observed 95% CL upper limit on S together with the parameterization of S as a function 
of TjG and 1 / Mg are shown in the left pane of Fig. [s] The intersection of the limit on S with 
the curves determines the upper 95% CL upper limits on the parameters rjc and 1/Mg. As 
seen from these plots, the upper limits are equal to 0.070 TeV^^ for t/q and 0.078 TeV^^ for 
1/Mg (for Wed = 2). We further translate these limits, by means of Eq. 0, into lower bounds 
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Figure 3: Parameterization of S as a function of the strength of the ED effect, //g for all cases 
except HLZ tied = 2 (left) and as a function of 1/M| for the wed = 2 case (right). The parame- 
terization shown is a fit according to Eq. (|5]). The solid line is the 95% CL exclusion limits on S, 
and is matched to the corresponding limits on t]Q and 1/M|. 



Table 3: 95% CL limits on Ms (TeV), as a function of the convention and number of ED. A 
comparison of the limits with a truncation of the production cross section above > Ms is 
also shown. The two limits for the Hewett convention correspond to positive and negative 
interference effects. 





GRW 


Hewett 








HLZ 










Pos. Neg. 


nED = 2 


riED = 3 


nED = 


4 = 5 


"ED = 6 


"ED = 7 


Full 


1.94 


1.74 1.71 


1.89 


2.31 


1.94 


1.76 


1.63 


1.55 


Trunc. 


1.84 


1.60 1.50 


1.80 


2.23 


1.84 


1.63 


1.46 


1.31 



on the fundamental Planck scale for various numbers of extra dimensions, tiee), as shown in 
Table |3] This is calculated trivially for Med = 2 and for Wed > 2 by using Eq. The limits in 
convention [5| are identical to the HLZ limits for tied = 4; the limit for the Hewett convention 
with constructive interference is 1.74 TeV and is close to the HLZ limit for Med = 5. 

We note that the LO signal cross section calculations become non-perturbative when the value 
of s in the 2 — > 2 process exceeds Mg. This effect is not taken into account in the SHERPA 
cross section calculations used in this analysis, or in previous studies of this process at the 
Tevatron [lOJ, where the effect is not expected to be important due to the lower collider energy. 
Because the energy of the LHC is significantly higher than the limits on Ms we are able to 
set in this analysis, we take this effect into account by conservatively assuming that the signal 
cross section is zero for \/§ > Ms- Under these assumptions, the limits on Ms decrease by 5% 
for nED = 2 (1.80 TeV) and 15% for Med = 7 (1.31 TeV). A summary of the limits under the 
assumption of a truncated production cross section is also shown in Table |3] 

In addition to setting limits on a specific model of large extra dimensions, we can also set a 
model-independent limit on any new physics model which results in central, high-£x dipho- 
tons, either resonant or non-resonant (e.g., Kaluza-Klein gravitons in the Randall-Sundrum 
model 1 25 1). We measure a 95% CL exclusion on the cross section times branching fraction 
times acceptance of 0.110 pb, for diphoton pairs with M^y > 500 GeV and the following kine- 
matic requirements on each of the two photons: Ej > 30 GeV and |?/| < 1.44. 
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6 Conclusions 



6 Conclusions 

In conclusion, we have performed a search for large extra dimensions ||TJ |2I in the diphoton 
final state with a data sample collected in pp collisions at v^s = 7 TeV corresponding to an 
integrated luminosity of 36 pb^^. We optimize the signal selection to reach maximum sensi- 
tivity in a counting experiment in a one-sided mass window by selecting events with centrally 
produced photons and large diphoton invariant mass. Given the absence of an excess over the 
SM direct diphoton background, we set lower limits on the cutoff scale Ms in the range 1.6-2.3 
TeV. While this analysis was being finalized, a phenomenological interpretation of the dijet an- 
gular distribution results from the CMS and ATLAS experiments appeared [26J and suggested 
even stronger limits on Ms. However, a dedicated experimental analysis and interpretation of 
the dijet data in the models with large extra dimensions has yet to be conducted. The results 
presented in this Letter extend the current limits reached at the Tevatron Q [lOl in all but the 
Med = 2 case. 
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